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Abstract. The cubature on Wiener space method, a high-order weak ap- 
proximation scheme, is established for SPDEs in the case of unbounded char- 
acteristics and unbounded payoffs. We first introduce a recently described 
flexible functional analytic framework, so called weighted spaces, where Fcller- 
like properties hold. A refined analysis of vector fields on weighted spaces 
then yields optimal convergence rates of cubature methods for stochastic par- 
tial differential equations of Da Prato-Zabczyk type. The ubiquitous stability 
for the local approximation operator within the functional analytic setting is 
proved for SPDEs, however, in the infinite dimensional case we need a newly 
introduced assumption on weak symmetry of the cubature formula. In finite 
dimensions, we use the UFG condition to obtain optimal rates of convergence 
on non-uniform meshes for nonsmooth payoffs with exponential growth. 



1. Introduction 

Cubature on Wiener space, a realization of the abstract KLV high order method 
after Shigeo Kusuoka [21], Terry Lyons and Nicolas Victoir [25] . is a weak ap- 
proximation scheme for stochastic differential equations. Significant advantages 
in comparison to other weak approximation schemes such as Taylor methods, see 
[TO] , are that it respects the geometry of the problem, and that at least theoretic- 
ally, it is possible to reach arbitrarily high rates of convergence without requiring 
the calculation of higher derivatives, see [25] Theorem 2.4, Proposition 2.5]. The 
concrete construction of such cubature paths of high order is still quite difficult, 
see |16) for paths up to order 11 for a single driving Brownian motion. Cubature 
schemes provide a time-discretization approximating the unknown expected value 
of a functional of the solution process of the SPDE by an expectation of an iter- 
atively constructed function on a high-dimensional discrete product space. Often 
a direct evaluation of the functional on the discrete probability space is too ex- 
pensive, therefore, several methods to speed-up the evaluation of cubature schemes 
such as recombination [53] GO] or tree-based branching [B] have emerged. Other- 
wise the functionals have to be evaluated with Monte Carlo or Quasi Monte Carlo 
algorithms on the discrete product probability space. 

High-order weak approximation schemes provide interesting lower complexity 
alternatives to standard multi-level Monte Carlo schemes if Quasi Monte Carlo 
algorithms or deterministic algorithms can be applied for the evaluation of the con- 
structed functionals. Indeed, multi-level Monte Carlo schemes lead to complexity 
estimates of order (almost) 0(e~ 2 ), i.e. to reach accuracy e a number of opera- 
tions of order e -2 is necessary. In contrast, QMC evaluations of weak, high-order 
approximation schemes of order k lead to complexity estimates of order (almost) 
0{e^ 1 ^ 1 / k ) 1 as long as the QMC integration yields optimal convergence (this in 
turn also depends on the dimension of integration space, which is moderate for 
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high order methods). Hence we believe that it is worth analyzing in depth the 
functional analytic framework of cubature schemes, i.e. we aim for constructing a 
flexible enough pool of Banach spaces of payoffs and Banach spaces of character- 
istics, where relevant problems from practice can be embedded. 

In this work, we shall relax the regularity assumptions of the cubature method, 
similarly as was done in [T2J [TOJ [T3] for the splitting approach of Syoiti Ninomiya 
and Nicolas Victoir [37]. Consider a stochastic differential equation on R™ in its 
Stratonovich form, 



All initial work was based on the fundamental assumption that the vector fields 
Vj : R n — ► R™ are bounded and C°°-bounded. This is also a typical assumption 
in other approximation methods for stochastic differential equations, e.g., in |33j . 
Some success in relaxing these assumptions, which are actually rarely satisfied in 
practical problems was achieved, at least for approximations of the splitting type, in 
works by Tanaka and Kohatsu-Higa [53] and Alfonsi pQ. While in the first one the 
focus was on extensions to Levy driving noise and in the second to CIR processes, 
in both approaches it was recognized that polynomially bounded payoffs are the 
correct context for problems with Lipschitz continuous vector fields. 

Another approach was suggested in [12]. There, splitting schemes were ana- 
lyzed on general weighted spaces, allowing in particular the approximation of Da 
Prato-Zabczyk stochastic partial differential equations where the drift part, the in- 
finitesimal generator of a strongly continuous semigroup on the infinite dimensional 
state space, is not even continuous. 

All these approaches profited from the special structure of splitting schemes, as 
there the stability or power boundedness of the discrete approximation operator 
can be shown by investigating every part separately. Instead, we follow a similar 
idea as was applied to the stochastic Navier-Stokes equations in [TT]. We extend 
the results of [3J to more general coefficients and payoffs. This allows us to obtain 
methods of order higher than 2 without having to resort to extrapolation, see [H [35] . 
While the use of the weighted spaces from [301 E2] is also mandatory here, we shall 
provide a refined analysis of the vector fields defined on these spaces. This will 
allow us to do a Taylor expansion of the cubature approximations to compute the 
local approximation order. 

While dealing with the stability, we shall use two different approaches. In the 
finite dimensional case with sufficiently smooth vector fields, the Gronwall inequal- 
ity yields the claim in a straightforward manner under a reasonable assumption 
of compatibility between the vector fields and the weight function. In the infinite 
dimensional case, we apply the method of the moving frame from [35) . This leads 
to time dependent vector fields that are nonsmooth in the time component. As this 
makes a Taylor expansion impossible, we introduce a weak symmetry condition on 
cubature paths, an assumption usually satisfied by cubature schemes. This allows 
us to obtain stability not only for Da Prato-Zabczyk equations with pseudocontract- 
ive generator, but also for stochastic differential equations on infinite dimensional 
state spaces, where the vector fields depend roughly, i.e., continuously, but not 
differentiably, on time. 

Finally, we consider the effects of the UFG condition in our setting. Under the 
same assumptions on the coefficients as in [5], we are able to prove optimal rates 
of convergence on non-uniform meshes for nonsmooth payoffs that are allowed to 
grow exponentially. 
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There are many successful discretisation schemes for stochastic partial differ- 
ential equations. |17j gives an overview of strong and pathwise schemes. Weak 
approximation schemes are more difficult. Recently, it was proved in [9] that an 
implicit Euler scheme converges almost with weak rate 1/2 for equations driven by 
space-time white noise, doubling the corresponding strong rate of convergence; see 
also the references in [pj for more background on weak approximation schemes for 
stochastic partial differential equations with space-time white noise. In contrast, 
we restrict ourselves to finite-dimensional driving noise, but obtain the same weak 
rate of convergence as for finite-dimensional state spaces. 

In our proofs, C denotes a generic positive real constant that can change from 
line to line. 

2. SPACES 

We recall the following definition of spaces of functions with controlled growth, 
see also [301 QH CLU QH EH • Notice that we obtain Feller-like properties for SPDEs 
in this setting. 

Definition 2.1. Let (X, ||-||x) be the dual space of a separable Banach space, 
and ip: X — > (0,oo) be bounded from below by some S > 0. For a Banach space 
(Y,|H|y), we set 

(2.1) B*(X;Y):=\f:X^Y: sup ^(a:)" 1 ||/(z)||y < oo) , 

I x£X J 

endowed with the (/3-norm 

(2.2) \\f\\ v := sup^ri/OrOHr. 

xex 

Let k > 0. If ip = (fj)j=o t ... t ki Pj '■ X (0, oo) bounded from below by some 5 > 0, 
j = 0, . . . , fc, we set 

B%(X;Y) := {/ G C k (X;Y): sup ^(x)- 1 ]]!^ f(x)\\ Lj[x . Y) < oo 

xex 

(2.3) for j = 0,...,k}. 
B"£(X; Y) is endowed with the norm 

h 

(2-4) ll/IU-ll/IL + El/Ui, 

where the seminorms | • | v . j are given by 

(2-5) := sup VjW^W&fWWLMY). 

xex 

Here, Lj(X; Y) denotes the space of bounded multilinear forms a: X^ — > Y, and is 
endowed with the norm 

(2.6) \\a\\ Lj {x-Y) ■= sup ||o(/li,...,/lj-)||y. 

||fci||<i,i=i,...,j 

For simplicity, we set Lo(X;Y) :— Y; we remark that L±(X;Y) is the space of 
bounded linear operators X — > Y, and in this case, the above norm is the usual 
operator norm. If Y = M, we define B^(X) := B^(X;R) and B%(X) := B%(X;R). 

Definition 2.2. Let (X, \\-\\x) be the dual space of a separable Banach space. A 
function ip is called admissible weight function if and only if tp : X —¥ (0, oo) is such 
that Kr := {x G X : (p(x) < R} is weak-* compact for all R > 0. 

It is called D- admissible weight junction if and only if it is an admissible weight 
function and for every x G X, there exists some R > such that B e (x) C Kr for 
some e > 0, where B e (x) := {y G X : \\y — x\\x £ e} is the closed e-ball around x. 
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It is called C-admissible weight function if and only if ip is bounded from below, 
weak-* lower semicontinuous, and if for every x <E X, there exists some e > such 
that (p is bounded on B s (x). 

Remark 2.3. We do not require C-admissible weight functions to be admissible. 
However, ip is D-admissible if and only if it is admissible and C-admissible. 

Theorem 2.4. Let k G N, and assume that ip = (tpj)j=o,...,k is a vector of C- 
admissible weight functions. Then, ~B^(X;Y) is a Banach space. 

Proof. Let (f n )n&i be a Cauchy sequence in this space. It is clear that /„ admits 
a pointwise limit /. Moreover, it follows that for every x € X and every closed 
e-ball B e (x), f n \B e (x) are Cauchy sequences in C k (B e (x); Y). But this entails that 
J\b c (x) € G k (B e (x); Y). As differentiability is a local property, we see that / e 
C k (X; Y). The necessary estimates for / and its derivatives are now easy to see. □ 

Remark 2.5. A counterexample showing the necessity of C-admissibility in The- 
orem 12.41 is given in Appendix [A] 

Definition 2.6. Let (X, \\-\\x) be the dual space of a separable Banach space, its 
predual being W, X — W* , and (Y, ||-||y) a Banach space. The space of bounded 
smooth cylindrical functions is defined by 

A(X,Y) := {f:X^Y:f = g((; Wl ),...,(;w n )) 

for some g £ G^(R n ;Y), 

(2.7) Wi e W, i = l,...,n, n € N}. 

Here, (•,•) denotes the dual pairing of X and W. For Y — R, we set A(X) :— 
A(X,M). 

Definition 2.7. Let (X, \\-\\x) be the dual space of a separable Banach space and 
(y, \\-\\y) be a Banach space. Let ip be an admissible weight function on X. 

The space B^{X\ Y) is the closure of A(X, Y) in B lp (X; Y). For Y = R, we set 
B*{X) := B^(X;R). 

Remark 2.8. [HI Theorem 4.2] shows that our definition of B^(X) here agrees 
with our earlier definition from [121 Definition 2.2]. Due to Theorem 2.7], the 
functions in B^(X) are characterized by the property that both /\k r € C((Kr) w *) 
and 

(2.8) lim sup ip(x)- 1 \f(x)\=Q. 

Definition 2.9. Let (X, \\-\\x) be the dual space of a separable Banach space and 
(Y, ||-||y) be a Banach space. Let ip = (ipj)j=o,...,k with ipj D-admissible weight func- 
tions for j = 0, ... , k. The space BfiX; Y) is the closure of A{X, Z) in Bf(X; Y). 
For y = R, we set := B%(X;R). In particular, by Theorem^ it follows 

that Bf(X) is a separable Banach space. 

One essential property of B^(X) spaces is that the dual space of this separable 
Banach space is a well understood space of Radon measures, such as in the case of 
Cq(X) for locally compact spaces X. 

Theorem 2.10 (Riesz representation for B^{X)). Let I: B^{X) -> R be a con- 
tinuous linear functional. Then, there exists a finite signed Radon measure fi on X 
such that 

(2.9) 1(f) - / f(xMdx) for all f E B*(X). 

Jx 
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Furthermore, 

(2-10) / t[>{x)\n\(dx) = \\e\\ L (B*(x),R), 

where denotes the total variation measure of \x. 

As every such measure defines a continuous linear functional on B^(X), this 
completely characterizes the dual space of B^{X). 

This allows for the introduction of the generalized Feller property, such that we 
can speak about strongly continuous semigroups on spaces of functions with growth 
controlled by ip, in particular functions which are in general unbounded. 

Let (P*)t>o be a family of bounded linear operators P t : B^(X) — > B^(X) with 
the following properties: 

(Fl) P = I, the identity on B*{X), 

(F2) P t+S = P t P s for all t, s > 0, 

(F3) for all / <= B*(X) and x £ X, lim t ^ + Ptf(x) = f(x), 

(F4) there exist a constant C £ K and e > such that for all t £ [0, e], 

||-P*IU(B*(Jf)) < Cj 

(F5) P t is positive for all t > 0, that is, for / £ B^{X), f > 0, we have P t f > 0. 
Alluding to jTSJ Chapter 17], such a family of operators will be called a generalized 
Feller semigroup. 

We shall now prove that semigroups satisfying IF1I to IF4I are actually strongly 
continuous, a direct consequence of Lebesgue's dominated convergence theorem 
with respect to the measure existing due to Riesz representation. 

Theorem 2.11. Let (Pt)t>o satisfv WT\ to lF4l Then, (Pt)t>o is strongly continuous 
on B^(X), that is, 

(2.11) t lim \\P t f-fU = for all f e B*(X). 

Proof. By [TU Theorem 1.5.8], we only have to prove that t i-> l(Ptf) is right 
continuous at zero for every / £ B^(X) and every continuous linear functional 
I: B^{X) — > M.. Due to Theorem l2.101 we know that there exists a signed measure 
i/onl such that £(g) — J x gdv for every g £ B^'(X). By |F41 we see that for every 
t£ [0,e], 

(2.12) \Ptf(x)\ < C1>(x). 

Due to (12.101) , the dominated convergence theorem yields 

(2.13) lim / P t f(x)v(dx) = [ f(x)v(dx), 
t-^Q+Jx Jx 

and the claim follows. Here, the integrability of ip with respect to the total variation 
measure |f| enters in an essential way. □ 

3. Vector fields and directional derivatives 

When we ask for convergence rates we have to specify large enough sets of test 
functions within the basic B^(X)-spaces. For this purpose we need to analyze 
directional derivatives and their functional analytic behavior. This can be done 
within the setting of B^(X; Y) spaces. 

Let (X, \\-\\x) be the dual space of a separable Banach space. Given (Z, \\-\\z) 
the dual space of another separable Banach space that is embedded in X, we derive 

conditions onV: Z — > X such that the directional derivative g £ Bt^Z), where 
(3.1) g(z) := Df(z)(V(z)) for z £ Z 
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and / € Bt (X). Here, ip = (ipj)j=o,...,k and tp = {^j)j=o,...,k— l are vectors of 
D-admissible weight functions on X and Z, respectively. 

We shall assume that V £ B^_ 1 (Z;X) for some vector <p = (<Pj)j=o,...,k-i of 
C-admissible weight functions on Z. Then, V is k — 1 times continuously Frechet 
differentiable. As / G C fc (X), the Leibniz rule yields 
(3.2) 

J 1 

D 3 g(z)[hi,-' ■ ,hj) = 2J ■ _ ^ X! 9j,ii z > h <nr- ■ ,h aj ), 3 = 0,...,k- 1. 

i=o o-eSj 

Here, Sj denotes the symmetric group with j elements, and 

(3.3) h u ■ ■ ■ , hj) := D t+1 f{z){h u ...,hi, Dt-'VWihi+i, • • • , fy)). 
In particular, if we assume that for some constant C > 0, 

(3.4) ^(z) >C- X W^Wi^Vj-^) for j = 0,— ,fc-l, 

i=o w 

it follows that 3 € V,\_ X {Z). 

It is not so straightforward to prove that g can also be approximated by functions 
in A(X), which would imply g £ Bf_ 1 (Z). In [TU], a general theory for multiplic- 
ation operators on spaces is derived. Here, we take a different route, focusing 
on the problem at hand. The following definition is essential. 

Definition 3.1. Given a Banach space (X, \\-\\x) and the dual space (Z, \\-\\z) of a 
separable Banach space. Let V £ B"£(Z;X) with ip a given vector of C-admissible 
weight functions on Z. We say that V £ C^(Z; X) if and only if for every y £ X* . 
there exists a constant Cy, y > such that for all R > 0, there exists a sequence 
(fn)neN C A(Z) with sup„ eN 1 1 v n \ \ v , k < Cy.y such that, with v := y o V, 

(3.5) n^L^ ~ Vn Wc k (B R (0)) = 0. 

Here, -Br(O) is the closed unit ball of radius R in Z, and 

fe 

(3-6) IIsIIc*(br(o)) : = Y] su p W&gWWLjW- 

Remark 3.2. It is clear that vector fields such as those from [T31 Section 2.2] satisfy 
the above assumption. More generally, if Z is a Hilbert space and is compactly em- 
bedded into a larger Hilbert space Y such that y o V can be extended to a smooth 
mapping Y — > K lying in C^(Y;M.) for all y £ X*, then the above assumption is 
satisfied, i.e., V £ C%(Z;X) for every vector tp of C-admissible weight functions 
on Z. Indeed, the extension of y o V and its derivatives are continuous on Y, 
whence uniformly continous on the compact set Br(0). Let us fix a sequence of 
increasing finite-dimensional, orthogonal projections 7r„ — > idy converging strongly 
to the identity: composing the extension of y o V with n n yields a pointwise con- 
verging, equicontinous sequence of cylindrical function on -Br(O), which is - up to 
a smoothing argument - the desired assertion. 

See also [THl Theorem 5] and [TQ1 Theorem 2.39] for comparable arguments. In 
particular, this implies that Nemytskii operators are included in our setup if Z is a 
Sobolev space of sufficiently smooth functions, see also [TO] Example 2.48]. 

This definition should also be compared to the form of the multiplicative noise 
suggested in [9l Remark 2.3]. It is similar in spirit to the definition of C£(H;H), as 
there, A is assumed to be a negative self-adjoint operator with a compact inverse. 
Hence, if we consider a single component of the noise, x i-> cr((— A)~ 1 / 4 x), with 
ct H — y H a C 3 -function with derivatives bounded up to order 3, it satisfies our 
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assumptions given above and hence lies in C^(H) with <po(x) := (1 + H^llff) 1 ' 2 and 
¥>j( x ) ■= l , 3 > 1- 

Theorem 3.3. Fix k > 1. Let ip = (V'i)i=o,....fe ^ e a vector of D- admissible 
weight functions on X, and ip — (ijjj)j=Q,...,k-i a vector of D- admissible and <p — 
{<Pj)j=o,...,k-i o, vector of C'-admissible weight functions on Z. Suppose (|3.4[) . 

Then, the Lie derivative C: Cu_i(Z;X) x Bt(X) — > Bt^Z) defined through 

(3.7) C(V, f)(z) := C v f(z) := Df(z)(V(z)) 
is a bilinear, bounded operator. 

Remark 3.4. Clearly, it is necessary that V € C^(Z; X) if Lyf € B^(Z) is supposed 
to hold for / 6 for a sufBciently large class of weight functions ip. Indeed, 

choose ipo(x) := /o(||x||x) with some increasing, left continuous and superlinear 
function p, and ipj arbitrary D-admissible weight functions on X. Then, /:=j/6 

Bf +1 (X) for all y € X*. Hence, £v/0) = y(F(z)), andi/o^e (Z) implies 
that 7eC[(Z;I). 

Proof. The claimed boundedness of £ was remarked above, and follows straight 
away from (|3.4p . 

Hence, we only need to prove that Cvf & ^t-i(^) ^ or gi ven ^ e ^k-i(^> X) and 
/ = toi), • • • , (•, w n )) £ -4(X); the result then follows from a density argument. 
Fix e > 0. We shall construct g e e A(Z) such that ||£y/ — g E \\^ k < Ce with some 
constant C > independent of e. 

Choose a dual set of vectors (Ci)i=i,...,n C 2 of K)i=i,...,nj i-e., (d,vjj) — Sij. 
Let Z„ := span {£,- : £ = 1, • • • , n}, and define ir : X — > Z n by irx :— ^ i=1 (s, Wi)Q. 
Then, / o 7T = /, and 

(3.8) £ v /(z) = ^I)/(*)(Ci)<nzW- 

i=l 

Clearly, to, € X* , and thus by Definition l3.ll there exists Cy := maxi = i.... in Cy. Wi > 
such that for all R > 0, we can find w}j e € -4(Z) with ||u]j e |L &_i < Cy and 

(3-9) |K o V - Vfl, e ||c»-i(-Ba(0)) < e > 

where £?i?(0) denotes the closed unit ball in Z. Setting g E :— J27=i Df(-)(Q)v l Re £ 
A(Z), it follows that with a constant C/ > independent of R > 0, 

(3-10) ||£v/-fl , e||c-i(B K (o)) < Cje. 

Choose R e > large enough such that ipj(z) > £ _1 for ||z||z > -Re- This is possible 
as the embedding Z — > X is continuous. Hence, as / and all its derivatives are 
bounded, 

(3.11) ^ j (z)- 1 \\D j C v f(z)\\L j ( Z )<C f e for \\z\\ z > R e , j = 0, . . . , k - 1, 
where C/ is independent of s. Furthermore, 

(3.12) ^(«)- 1 ||D i fl E («)|U,(z ) < C f , V E for \\z\\ z > R,, j = 0, . . . , k - 1, 

where Cfy > depends on / and V, but not on e or R e . Plugging the results 
together proves the claim. □ 

Let us consider two special cases. 
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Corollary 3.5. Let (H, \\-\\h) be a Hilbert space, (Z, \\-\\z) o, continuously embedded 
Hilbert space. Define the D-admissible weight functions ipj( x ) '■— C0Sn (IMI.ff) on H 
and tpj( x ) : ~ cosh(|| z) on Z and the C- admissible weight functions ipj(x) '■= 1 
on Z, j > 0. Then, for every k > 0, the mapping 

(3.13) L:Bt{X)xCt_ 1 {Z;X)^Bt 1 {Z), (f,V)»£ v f, 
given by Cy f(x) :— Df(x)V(x), is bounded and bilinear. 

Remark 3.6. If Z = H, this has the simple interpretation that bounded vector 
fields map cosh- weighted spaces into themselves. 

Proof. This is straightforward from Theorem 13.31 as the defined there is only a 
multiple of ipj in this case. □ 

The following special case is very useful in the analysis of stochastic partial 
differential equations of Da Prato-Zabczyk type. 

Corollary 3.7. Let (H, \\-\\h) be a Hilbert space, (Z, \\-\\z) o, continuously embedded 
Hilbert space. Fix n £ N. Define the D-admissible weight functions ipj( x ) := 
(1 + jWI^)^')/ 2 on H and ^(x) := (1 + |M||) ( "^ )/2 on Z, j = 0, . . . , n - 1, 
and the C-admissible weight functions (fo(x) := (1 + H^Hfj) 1 ' 2 and fj(x) := 1 on 
Z , j £ N. Then, for k < n — 1, the mapping 

(3.14) C:Bt(X)xq_ 1 (Z;X)^Bt 1 {Z), (f,V)^C v f, 
given by Cy f{x) :— Df(x)V(x), is bounded and bilinear. 

Remark 3.8. This means that linearly bounded vector fields Z — > X with bounded 
derivatives (hence also Lipschitz continuous) map polynomially bounded functions 
to polynomially bounded functions, with the same weights. In particular, if A : dom A C 
H — » H is a densely defined, closed operator, then Va £ Cj^ (dom A; H) for all k > 0, 
where (p is defined as in Corollary 13.71 dom A is endowed with the operator norm, 
and Va(x) := Ax for x £ dom A. 

Proof. Calculating 

(i + \\x\\%)w*(i + iwiD 1 / 2 + £ (0 (i + Mir- 1 -^ 2 

(3.15) <Ciij(x), 

the claim again follows from an application of Theorem 13.31 □ 

4. Stability of cubature schemes 

We shall now prove stability of cubature on Wiener space in the setting of 
weighted spaces. Consider from now on the following setup. Let on [0, 1] be 

given paths (u}f')i=i,...,N, (s) = (^i j (s))j=o,...,d, w| 1),0 (s) = s, and weights 
(^i)i=i,...,N of a cubature on Wiener space of order m > 1 for a <i-dimensional 
Brownian motion, i.e., for all multi-indices . . . ,jk) with k + # {i : ji = 0} < m 
and a d-dimensional Brownian motion (B\ )j=i,...,d,4>0j 



(4.1) E 



o&Bi\ ■ ■ • o dBH 

0<si<---<s fc <l 



N 

X A / • • / d^ 1)jl ( Sl )...od W p fc ( Sfe ). 
' J0<s 1 <-<s k <l 



i=l 
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Here, we have set := t and odB a := dt for ease of notation. For a general time 
interval [0, At], we set 

(4.2) 4 At) '°( s )~ s and 4 At) ' j ( s ) :=V^uP' S (a/At), j = l,...,d, 

so that {u>i AtS> )i=i,...,N and (Xi)%=i,...,N define a cubature formula on Wiener space 
of order m on [0, At]. The approximation of the Markov semigroup (Pt)t>o, given 
by P t f(x) := E[f(X?)] for a function f:H->M., where (X?) t > solves the Stra- 
tonovich stochastic differential equation 

d 

(4.3) dX? = '£v j (X?)odBi, X$=x, 

j=o 

on some state space H, then reads 

(4-4) P t f(x) « Qf t/n) f(x), 

where the one step approximation operator is defined by 

N 

(4.5) Q(A t )/(a:):=E A «/( X A t (^ At) )), 

»=i 

with Xf{uj^') the solution of the problem 

(4.6) dX s >| At >) = £ VjiX^^d^is), X >j At) ) = *. 

3=0 

Under certain smoothness assumptions on the vector fields Vj, j = 0, . . . , d, and 
the payoff /, we expect that 

(4-7) \P f f(x) - Q n [t/n) f{x)\ < Cn-^' 2 , 

where the constant C > can depend on /, Vj, j = 0, . . . , d, and x € H. For the 
case H finite-dimensional and / and Vj bounded and C°°-bounded, j = 0, . . . , d, 
it is known that C depends on the supremum norms of / and its derivatives, but 
not on x € H, see [3S]. For more background on the method, see [23 13 |3]. An 
alternative approach can be found in [21] [22]. Its implementation as a splitting 
method is given in [27] , see also [3S] [T] [33] . 

Our strategy is as follows. First, we consider the finite dimensional case. Here, 
the analysis is straightforward. Afterwards, we turn to the infinite dimensional 
setting. Here, our aim is to prove stability for Da Prato-Zabczyk equations with 
pseudodissipative generator. We prove first the auxiliary result in Theorem 14.41 
which might be of independent interest. The method of the moving frame then 
yields first Theorem 14 .71 and the Szokefalvi-Nagy theorem allows us to conclude in 
Corollary EE] 

4.1. Finite dimensional state space. Given a Stratonovich SDE on K", 

d 

(4.8) dX? = Y,Vj(X?)odBt, X$=x, 

3=0 

we let the local discretisation of P t f(x) := E[/(Af )] be defined by 

N 

(4-9) Q (At) /(x):=EAJ(A At ( Wi (At) )), 

i=l 



10 PHILIPP DOR.SEK, JOSEF TEICHMANN, AND DEJAN VELUSCEK 

where X^{lo\^) is the solution of the problem 

(4.10) AX*M M) ) = !>(*.>J At) )W At)j '(*)> *oK (At) ) = 

Theorem 4.1. Let tp be an admissible weight function on K n , and assume that 

(4.11) \ViVjip(x)\ + \Viip(x)\ < Cij){x) fori = 0,---,dandj = l,---,d, 

where we require that all the necessary derivatives are well-defined. 
Then, there exists a constant C > independent of At > such that 

(4.12) Q(At)V»(aO < exp(CAt)ij(x). 
Proof. We define the intermediate operator 

JV 



(4.13) Q(At,s)m '■= E A ^W(^ (At) )) for s E [0,t] 

i=l 

and note that Q(At) = Q(At,At)- The definition of the iteration step yields 

(4.14) i,(X*(4 At) ))=^(x)+J2 f S VM^(4 M) ))d4 At) ' j (r) 



^(x) + f yo^(^K (At) ))dr + ^ ^xH (At)j ( S ) 



j=l fc=0 

By em), 



, =lfc=0 ^o Jo 



(4.15) / l/ ^(^K (At) ))dr < C / ^(X r >f 4) ))dr 



Furthermore, as \u[ At),i {s)\ < C(At) 1 / 2 and |#^ At),;, '(s)| < C(Ai)- 1 / 2 , Fubini's 



^ (At),j / \ 



dg 



theorem yields 

^y^(^(^ At )))d^ (A4 ^(g)d^ At)j (r) 

<c/Vf*)-( S )-^-(g)|^(^(o;| At ))) 
Jo 

(4-16) <C f'i>{X^ At) ))dq. 

Jo 

Thus, we see that 

N 

(4.17) Q (A i,.)V(aO = E A 4 ^(^K (At) )) 

i=l 

< + E E A ^i At) ' J ( S ) + C / 0(At,r)^(x)dr. 

Defining «aj.s(j;) := Sj=i ^j'VK 3 ') Si=i ( s )> ^ ne Gronwall inequality yields 

(4.18) Q(At,s)^{x) < ip(x) + a At<s (x) + [ (^(x) +a Atir (x))Cevp(C(s-r))dr. 
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Note that a A t,At(x) = by the equality J2iLi Ai^f A *' J (At) = 0. Furthermore, 

C C 

(4.19) a AttS (x) < CVKiip(x) < — (1 + M)i/)(x) < — exp(At)tp(x). 

This proves 

Q(At)^{x) = Q (AtAt) ^(x) < i){x) (\ + ^exp(Ai)(exp(CAi) - l)j 

(4.20) < exp{CAt)ip(x), 

where C :— max(C 2 /2, C + 1), which is the required estimate. □ 

4.2. Time-dependent stochastic ordinary differential equations on Hil- 
bert space. Let H be a Hilbert space, and consider the nonautonomous stochastic 
ordinary differential equation 

d 

(4.21) dXf = £^(t,*?)odfl£, Xg=x, 

3=0 

on H . We define cubature approximations of (14.21(1 by 

(4.22) dX^( W ^) = X;^(* + *,^ ,x (^ At) ))dwi At)j \«), *r =*> 
and the approximation operator by 

iV 

(4.23) QUt)f(x) ■=T, X ^ X M(4 At) ))- 

i=l 

Definition 4.2. A cubature formula (w^ , Ai)i = i,.. jv is called symmetric if for 
every z g {1, . . . , AT}, there exists some i' € {1, . . . , TV} such that Ai = A-j/ and 

(4.24) ^ At) ' J {s) = -4 A *' J ( s ) for a11 s e [0, At] and j = 1, ... , d. 
It is called weakly symmetric if for j = 1, . . . , d, 

N 

(4.25) X M At) ' J ( s ) = fOT s e [0, At]. 

i=l 

Remark 4.3. Clearly, all symmetric cubature formulas are also weakly symmetric. 
Note that many known cubature formulas are actually symmetric. Moreover, a 
non-symmetric cubature formula can be made symmetric by adding the negatives 
of the paths with the same weights to it and dividing all weights by two. This will 
at most double the number of paths. Thus, if we use a cubature formula with a 
small number of paths in high dimensions, we can also find a symmetric cubature 
formula with this property. 

Theorem 4.4. Suppose that the cubature formula used in the definition of Q* A t) 
is weakly symmetric. Let tp be an admissible weight function on H and suppose 

(4.26) \\Dij{x)\\<C{\ + \\x\\ 2 )- 1/2 ^{x) and 

(4.27) \\D 2 ij{x)\\ < C(l + Wxfy^ix) 

with some constant C > 0, Furthermore, assume that for some constant C > 
independent of t, 

(4.28) \\V.j{t,x)\\ < C(l + ||x|| 2 ) 1/2 forj = 0,...,d,xeXandt£[0,T], 

and that x i— > Vj(t,x) is continuously differentiable with derivative bounded uni- 
formly intE [0, T] for j = 1, . . . , d. 
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Then, there exists a constant C > such that for all t £ [0, T] and At £ [0, T— t], 

(4.29) Q\At)^( x ) < exp(Ct)i;(x) for all x <= H. 

Remark 4.5. The above result is remarkable as we do not assume that the vector 
fields Vj are differentiable with respect to t. This is also the reason why we cannot 
simply apply Theorem 14. II to conclude. 



Proof. Define the intermediate approximation for s £ [0, At] by 

N 

(4-30) QtAt, s) m ■■= E Ai/(X s *.>f 

i=l 

As in the proof of Theorem 14. 11 we note that Q' At ^.t) = Q^At)- For — s — At, 

*/«>( At) )) = 

, ,(At)» ,(At),j, 



(4.31) V(z)+E / ^(^l^l^^aJl^^d^^r), 

Consider <?j(r", x) := Dip(x)Vj(t + r, x). Then, 
g J (p,X^(J l At) ))=g J (p,x) 

(4-32) +E f X*->J At >))V k (i + g, ^( W | A {) ))do;f 



From (|4T26|) . (|4~27| and (|4T28|) . we obtain that for < s < At < T, 

\g {r,x)\ = \Dip(x)V (t + r,x)\ < C\\D^{x)\\ ■ ||V (t + r,a:)|| 

(4.33) < Ctp(x). 

We argue in a similar manner for D x gj(r, x)Vk(t + q,x), j = 1, . . . , d, k = 0, . . . , d, 
to obtain that for < q < r < At, 

\D x g 3 {r, x)V k (t + q,x)\ = \D 2 ^j{x){V 3 {t + r, x),V k (t + q,x)) 

+ Dil>(x)D x Vj (t + r, x)V k (t + q,x)\ 

(4.34) < Ctp(x). 

An application of Fubini's theorem just as in the proof of Theorem I4.ll gives 

^Ptf>j At) )) = m + f g (r, X^(4 At) ))dr + E f 9 S {r, x)d^ (At) <» 

Jo j =1 Jo 

+ EE f f D ^(r, X?{ U M))V k {t + q, X^^))M^%)M^Hr) 
j= ik=o Jo Jo 

(4.35) 



< 



^{x) + C ^K (At) ))dr + E / gj^x^^ir), 
Jo j=1 Jo 

where we apply that At < T. As from the weak symmetry of the cubature paths, 

N d fS d fS / N 

$>£/ ft(r ! ^ At )^(r)=E/ ftMKEW 4i,J w 

(4.36) = 0, 
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we obtain 

(4.37) Q{At,s)i>(x) < %l>{x) + C f Q(At,r)^(x)dr. 

Jo 

An application of Gronwall's lemma yields Q(At)i } { x ) < exp(CAt)i/j(x), which 
proves the result. □ 

Remark 4.6. It is clear that the given assumptions on the vector fields and the 
weight function are not the only ones possible. Instead, we could also require the 
vector fields to be bounded uniformly in t € [0, T], and allow the weight function 
to satisfy ||£)^(x)|| + ||D 2 '(/'(a;)|| < Cij){x). While the situation above corresponds 
to polynomially growing weight functions and linearly bounded vector fields, this 
variant corresponds to exponentially growing weight functions and bounded vector 
fields, see also Corollaries 13.71 and 13.51 

Such an approach might be more appropriate when dealing with exponentials of 
stochastic processes such as Levy processes, which are ubiquitous in applications in 
mathematical finance, as they ensure nonnegativity in a simple manner and allow 
us to work on the natural scale of the problem. 

4.3. Da Prato-Zabczyk equations. Suppose now that 

d 

(4.38) dXf = AX?dt + Y,V j (X?)odBi, X% = x, 

j=o 

is a stochastic partial differential equation of Da Prato-Zabczyk type on some Hu- 
bert space H, see [71 [5] for a comprehensive exposition of the theory of such equa- 
tions. Here, solutions are understood in the mild sense, 



Xf = exp{tA)x + ]T f . 

„-_n JO 



(4.39) Xf = exp{tA)x + ^ / exp((* - s)A)V j {Xf) o dB>, 

j=o- 

and we also define the cubature discretisations in the mild sense, 

(4.40) X?(J t At) )=eMtA)x + J2 f «p((t - s)A)V j (Xf 4) ))da;f ^(s). 
Again, the approximation of the Markov semigroup P t f(x) := E[/(Xf )] is given by 



N 



(4.41) Q(A*)/(z) := V Kf{XU4 At) ))- 



i=l 



Theorem 4.7. Suppose that A is the generator of a group St — exp(tA), t G M, 
and that the cubature formula used in the definition o/Q^t) is weakly symmetric. 
Let %j) be an admissible weight function on H. With some constant C > 0, let 
tjj(Stx) < exp(Ct)tj)(x) for all x € H and t > 0, and 

(4.42a) 11^^(^)11 < C(l + ||a;|| 2 ) _1/2 -0(^) and 

(4.42b) \\D 2 ip{x)\\ < C(l + ||a;|| 2 )-VW- 

Furthermore, assume that 

(4.43) \\Vj(x)\\ ^Cil + WxW 2 ) 1 / 2 forj = 0,...,d, 

and that Vj is continuously differ entiable with bounded derivative for j = 1, . . . , d. 

Then, for any T > 0, there exists a constant C > such that for every At € 
[0,T], the operator Q (At) satisfies 

(4.44) Q(At)i>(x) < exp(CAt)^(x) for all x G H. 
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Proof. We apply the method of the moving frame from [35]. This yields that 
Xf — S t Y t x , where (Y t x ) t >o satisfies the Hilbert space stochastic ordinary differen- 
tial equation 

d 

(4.45) dl? = X;^(*,57)°d^, y V = V, 

3=0 

with Vj(t, y) — S- t Vj(S t y)- Thus, rewriting the cubature discretisations of {Xf ) t >o 
using {Y t x )t> , 

(4.46) dF s >| At) ) =J2Vj&Ys^ M) ))M M),j (^ 

3=0 

we see that, if we define 

N 

(4-47) Q(A t )/(2/):-E A ^( r A t K (At) )) 

i=l 

for /: H — » M, then Q(At)h{x) = Q(At)9{x), where g{y) := h{SAty)- I n particular, 

(4.48) Q(M)1>(x) = Q(At)(^ ° < exp(CAt)Q (At) ^(a;), 

where we apply the assumptions on and the positivity of Q(At)- 

But now, we are in the situation of Theorem 14.41 The estimates for ip are clear 
by assumption, and for Vj{s,y), we note that, as s € [0,T], 

(4.49) ||V-( S , y)\\ = WS-VjiSsy)]] < C{1 + M 2 ) 1 ' 2 for j = 0, . . . ,d 
and 

(4.50) \\D y Vj{s,y)\\ = \\S- a D y V J (S B y)S a \\ < C for j = l,...,d. 
An appeal to Theorem 14.41 yields 

(4.51) Q(AtMx) < exp(CAt)^(x), 

and the result follows. □ 

The Szokefalvi-Nagy theorem now allows us to obtain a corresponding result for 
pseudocontractive semigroups. 

Corollary 4.8. Suppose that A is the generator of a semigroup of pseudocontrac- 
tions St = exp(iA), t > 0. Let ip{x) = p{\\x\\ 2 ) with some increasing and left 
continuous function p: [0, oo) — > (0,oo) (see also |12[ Example 4.1]J which satisfies 
p{Cu) < C p{u) for all u > and C > 0, and which is twice differentiate and 
satisfies 

(4.52) p'{u) < C(l + u)"V(w) and p"{u) < (7(1 + u)- 2 p{u). 

Furthermore, assume that ||^(.t)|| < C{\ + M 2 ) 1 / 2 for j = 0,...,d, and that Vj 
is continuously differentiable with bounded derivative for j = 1, . . . , d. 

Then, for any T > 0, there exists a constant C > such that for every At € 
[0, T], the operator Q (At) satisfies 

(4.53) Q{At)^{x) < exp{CAt)*P{x) for all xeH. 

Proof. Assume without loss of generality that {St)t>o is a semigroup of contractions. 
By the Szokefalvi-Nagy theorem [29l p. 452, Theoreme IV] , we see that we can find a 
Hilbert space {H, \\-\\-h) containing H as a closed subspace and a strongly continuous 
group {St)tem of unitary mappings such that St = irSt, where ir: H — > H is the 
orthogonal projection. 
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Define Tpuiu) '■— PdMlfJ an d V^iu) '■= Vj(ny), then it is easy to see that the 
assumptions of Theorem 14 .71 are satisfied. The results of [35] prove that the solution 
of 

(4.54) =s ty + J2 f S t _ s Vf (X s ^) o dBi 

satisfies Xf — nX t ' , and similarly for the cubature approximations. Setting 

N 

(4-55) QfaM E Ai/(X«f(4 At ))), 

i=i 

Theorem 14.71 yields that Q^^uiy) < exp(CAi)V>«(2/), and from i/)-h(%) = iA x ) 
for x € H we obtain that for x G H, 

N N 
i=l i=l 

(4.56) = Qf At) i>H{x) < exp(CAt)ijj n (x) = cx P (CAt)^(x). 

The result is thus proved. □ 

5. Convergence estimates of cubature schemes 

We are now ready to prove rates of convergence for cubature on Wiener space 
on weighted spaces. We shall only prove these results in the infinite-dimensional 
setting; corresponding results in finite dimensions are obtained in a similar manner. 

Let if be a Hilbert space and A the infinitesimal generator of a strongly con- 
tinuous semigroup of pseudocontractions on H. Fix £q E N. For £ — 0, . . . ,£o, let 
Hp be subspaces of H endowed with Hilbert norms H-Hffj, Hq = H, such that for 
£ = 0, . . . ,£o — 1, H( + i C H( and A: He+i —> He is a bounded linear operator. On 
He, we define D-admissible weight functions 

(5.1) := (i + iMiy s/2 , s>i, i = o,..., to, r--=r , 

and the functions 

(5.2) w t0 (x):=(l+\\xf H y /2 , w tj (x):=l, j>l. 

Define the vectors of weight functions ip^ :— (V'™ - J )j=o,...,fc> k < n, and <pt :— 

(<Pe,j)j=o,-,k- 

Assumption 5.1. The vector fields satisfy 

(5.3) Vj eC% e {H e ,H e ) for j = 0, . . . , d and £ = 0, ...,£ . 

Remark [XU shows that A G C^ +1 {H e+1 , Hi) for I = 0, . . . ,£ - 1. For z G i^, 
£ = 0, . . . , £oi we can then consider the Da Prato-Zabczyk equation 

d 

(5.4) dXf = AXfdt + ^^(Xf)od J Bf, X * = z, 

on As the assumptions on the vector fields Vj essentially mean that they 

are Lipschitz continuous with bounded derivatives, all these equations have unique 
solutions, agreeing if we vary £ for sufficiently smooth initial conditions. 

Assumption 5.2. The Markov semigroup (Pt)t>o, -Ft/ 0*0 := E[/(X t x )], is strongly 
continuous on B^ e (He) for all n G N and I = 0, . . . , Iq. For some fco G N, Pt is a 

bounded map from B k e (He) into itself for k = 0, . . . , ko and n G N, n > k, with 
norm bounded uniformly in t G [0, T] for every T > 0. 
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See also (T2J Section 5, Lemma 7.19] for sufficient conditions for these assump- 
tions. 

5.1. Taylor expansion of stochastic partial differential equations. 

Theorem 5.3. Let I — 1, ...,£q. Consider the strongly continuous semigroup 
(Pt)t>o on the space B^ e (Hi) with n>4. Denote its generator by (Q,domQ). 

Then, B 2 e (Hi-\) C domC?, and 

d 

(5.5) gf(x) = Df(x){Ax)+C v J{x) + -^^./W 

i=i 

for f <= BziHt-x) and x e H e . 

Proof. By the Ito formula [51 Theorem 7.2.1], it follows that for / 6 A(H(,-x), we 
have / G AomQ and / satisfies (|5.5[) . Corollary 13. 71 shows that the right hand side 

of (|5.5I) is a continuous linear operator 1 (H(^i) — ► (fli). The closedness 
of proves the claim. □ 

The next result follows directly from Corollary 13.71 together with the explicit 
representation in 



Corollary 5.4. Let k > 0. Under the assumptions of Theorem \5.S[ the infinites- 
imal generator Q satisfies the mapping property 

(5.6) Q:B^{Ht^)^Bt^{H t ) t £ = l,...,£ - 

Induction now yields: 

Corollary 5.5. Let j = £,..., £q. Under the assumptions of Theorem 1 5. 31 the 
powers of the infinitesimal generator Q satisfy 

(5-7) g>:Bt£ 3 (H e ^)^Bf\H e ). 

They are given explicitly by taking the powers of (|5.5[) . 

This allows us to obtain a Taylor expansion of Ptf for smooth enough /, which 
we will compare to the Taylor expansion of cubature approximations. 

Corollary 5.6. Let f e Baffc+x) ( H e-(k+i))> k + 1 < £ < £ 0) n > 2{k + 2). 
Then, 

k 

(5.8) Ptf = J2^ gj f + tk+lR ^ 

3=0 J ' 

where the linear operator Rt : k'- ^2(k+i) 1 ' > C^-(fc+x)) — * i^i) satisfies 

(5.9) \\Rt,kfh 7 <C T ||/||^_ Wiv2(fe+1) forte [0,T] 
/or a constant Ct > independent of f . 

5.2. Taylor expansion of cubature approximations. For a multiindex a = 
(ii, . . . , ifc), we define deg(a) := {j — 1, . . . , k: ij — 0}. The empty multiindex 
is denoted by 0, corresponds to k = 0, and satisfies deg(0) = 0. We set 

(5.10) A m := {a : deg(a) < m} and A* m := A m \ {0, (0)} . 
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Theorem 5.7. Assume that the cubature formula is of order m = 2k + 1. For 
f e B^l\f(H £ _ (k+1) ), k + l<£<l ,n>2(k + 2), 

(5.H) Q(At)f = E ^f"^/ + (Ai) fe+1 iW/, 

where the linear operator B,At,k ■ ^(fc+i)" 1 ' C^— (fc+i)) — * (-^) satisfies 

(5.12) < C T||/||^) fc+i)i2(fe+1) V At G [0,T] 
/or a constant Ct > independent of f. 

Proof. Under the assumptions on the vector fields, we can easily see that for every 
/ G w4(i?£_(/j + i)), we have the Taylor expansion 

(5.13) f(XlM M) )) 

= E ^ • • • ^/(*)4 1,,,,,< * ) («i (At) ) + #A tl */(*), 

(i I ,...,i*)e-4 ra 
where we define the iterated integrals by 

(5.14) lt-' ik) (^ At \g) 

0<ti<— <tk<At 

7 (ii,...,i»)^(At)j ;= f (ii,...,i fe )^(At) ; ^ the remainder term ^ t>fc / satisfies 

(5.15) R\ t<k m = E 4 0,, " ,<fc) (^ At) ,/( <0 ,...,i*)) > 

(tl,...,**)e^m 
(iO.il. ...,ifc)^ m 

and we set /3 (x) := Ax + V Q (x), /3j(x) := Vj(x), j = l,...,d, and /(i ,...,i fc ) := 
Pi . . . (3i k f , (io, G {0, . . . , <i} fe+1 . Summing up, it is easy to see by the 

scaling of the cubature paths that we can find a remainder term (At) k+1 as in the 
claim of the theorem with the correct estimates. To see that the initial terms have 
the form given, we use the order 2k + 1 of the cubature and the explicit formula of 
Q from Theorem 15.31 A density argument proves the result. □ 



5.3. The rate of convergence. We can now present our main result. 

Theorem 5.8 

2(k + 1) < fc 0j 



Theorem 5.8. For f G S^j 1 ' (H e _ (k+1) ), k + 1 < £ < £ 0) n > 2{k + 1), 



(5.16) \\P T f - Q? T/n) fU f < C T n-%f\\^ k+i)Mh+1) 

with a constant Ct independent of f. 

Proof. The local estimate follows from a combination of Corollary 15.61 and The- 
orem [3771 The stability of Q(T/ n ) horn Corollary 14.81 and the assumed invariance 

of ^(fc+i) (H£-(k+i)) with respect to P t prove the claim. □ 

Example 5.9. Let Hi — M. N be finite-dimensional, and assume n > 5; in the finite- 
dimensional setting, we do not need to consider subspaces of the state space. Then, 
/ G Bf n) (H) for all k > 0, where f(y) = y™, m = 1, . . . ,4. This implies that not 
only the expected value and the variance, but also the skewness and kurtosis are 
accurately computed by our scheme. Similarly, mixed moments are determined to 
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high accuracy, and if n is even larger, this also holds true for higher order moments. 
Such a property is very useful in risk management, where high precision in higher 
moments means an accurate evaluation of risk. Similar observations were made in 

CD HO- 

To illustrate this practically relevant behavior, we consider the Heston model, 
i.e., {X? , V^)t>o solves the stochastic differential equation 



(5.17) 
(5.18) 



ax? 

dV7 



k(6 - V t v )dt + (3y/V^d(pdBl 



v v 



and the stock price is given by S t := cxp(Xf), see, e.g., [2J. V t v is the squared 
stochastic volatility. As {Xf , V t v )t>o is a polynomial process, analytical formulas 
for the moments are available. For our simulation, we choose the parameters fi = 
.02, k = 5., 6 = .09, j3 = .6, and p = —.8. The starting values are chosen to 
be x = log(9.) and v = .0625, such that So = 9. We are interested in finding 
the mean m := E[Xf], the variance var := K[(X? — to) 2 ], the skewness skew := 
var~ 3 / 2 E[(Xf — to) 3 ], and the kurtosis kurt := var~ 2 E[(Jff — m) 4 ]. With the 
parameters given above, the exact values are found to be 



(5.19a) to = 2.192936688809, 

(5.19b) skew = -0.885007761283, and 



var = 0.019329503330, 
kurt = 4.321997672912. 



In our numerical simulation, we choose the cubature paths by applying a split- 
ting of Ninomiya-Victoir type, where the normal random variables are replaced by 
Gauss- Hermite quadrature, see also [34], where this splitting is also considered. For 
two-dimensional Brownian motion, this implies that we have 2q paths per interval, 
and that for i = 1, . . . , q, 



(5.20) 



(5.21) 



w S+i( s ) : = { 



(3s, 0,0 
(l,3&,i(s-l/3),0) J 
(1,6,1,3^,2(5-2/3)) 

(0,0,3&, 2 s) T , 
(0,3&,i(s -1/3), 6,2 
(3(s -2/3), ^,£,2) 



se [0,1/3], 

se [1/3,2/3], and 
se [2/3,1], 

[0,1/3], 

se [1/3,2/3], 
se [2/3,1]; 



see also [131 Example 2.2] for a similar rewriting of splitting-up methods. Here, 
(Cii Wi)i=i,...,q defines a g-point integration rule for a standard normal random 
variable in two dimensions of order 5. The weights of the cubature formula are 
Ai = \q+i = "Y". This cubature formula is weakly symmetric, as X)?=i w i£i,i — 
S?=i w i€i,2 = 0. The resulting ordinary differential equations can be solved ex- 
actly, see [231 [5] for the analytic formulas. 

For simplicity, we assume that (£,-, io,-),- = i r .. )g is the tensor product Gauss- Hermite 
quadrature of order 5. This implies q = 9, and hence, we have 18 paths per time 
interval. More efficient quadratures, in particular for the case of a high-dimensional 
driving Brownian motion, can be found in (32) . 

The numerical results are given in Figure [T] We observe the second order rate 
of convergence. With 8 time steps, we obtain a relative error of less than 10 -2 for 
all quantities of interest. 
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Figure 1 . Convergence of mean, variance, skewness and kurtosis 
of the log price in the Heston model 



Example 5.10. The Heath-Jarrow-Morton framework is included in our setup. As 
explained in [13] . it is more natural to use cosh- weighted spaces instead of poly- 
nomially weighted spaces in this case. The more general definition of vector fields 
in Definition 13.11 allows us to enlarge the class of admissible equations considerably 
compared to j!3j . 

5.4. Smoothing effects under the UFG condition. It is proved in [351 
that under the UFG condition, we obtain the optimal rate of convergence even for 
nonsmooth payoffs / by using non-uniform time steps due to the smoothing effects 
of Ptf in the direction of the vector fields Vj. The aim of this section is to show 
that a corresponding result holds true for unbounded payoffs, as well. In particular, 
we will focus on exponentially growing payoffs through the choice of cosh(a||a:||) as 
the weight function. This has important applications in mathematical finance, 
where one frequently models the log price as the solution of a stochastic differential 
equation, and thus, all payoffs will be a function of the exponential of the stochastic 
process. 

We assume that we are in the finite dimensional situation, H = l w for some 
N e N, and that A = 0. Suppose that all vector fields Vj : R N ^ R N are bounded 
and C°°-bounded. We choose the D-admissible weight function ip(x) := cosh(a|ja;|j ) 
for some a > 0. 

Proposition 5.11. For any a > 0, there exists C > such that 

(5.22) E[cosh(a||Xf||)] < exp(Ci) cosh(a||x||). 
Proof. Note that for any fceN, 

fc 

(5.23) L> fc cosh(a||a;||)(/ii,...,/ifc) < C k cosh(a||x||) JJll^ll- 

j=i 

With Q the generator of P t , we thus obtain from the boundedness of the vector 
fields that Q cosh(a||x||) < Ccosh(a||x||). It follows that 

(5.24) E[cosh(a||Xf||)] = cosh(a||x||) + f E[G cosh(a|| A7||)]ds 

^o 

(5.25) < cosh(a||x||) + / CE[cosh(a||X s x ||)]ds. 

^o 
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The Gronwall inequality now proves the claim. □ 

Corollary 5.12. For any p £ [1, oo) and T > 0, there exists Ct > such that 

(5.26) E[cosh(a||Xf||) p ] 1/p < C T cosh(a||x||) for all t £ [0, T]. 

Proof. We only need to note that for any p € [l,oo), there exists some constant 
C > with C^ 1 cosh(pu) < cosh(u) p < Ccosh(pu) for all u £ [0, oo), and apply 
Proposition [nHH □ 

We formulate now the ellipticity assumptions that are necessary to obtain smooth- 
ing effects. We follow [5]. 

The UFG condition: There exists I g N such that for every a £ A* , there 



exist tp at p e C^°(R Ar ), j3 € A* t , such that 



(5-27) V [a] = V>a,fiV w . 

PeA* e 

The VO condition: For some tpp £ C^{R N ), ft £ A%, 
(5.28) V = <PP V W' 

Theorem 5.13. Assume that the UFG and VO conditions are satisfied. Then, for 
any f £ C^°(IR Ar ), any k, m > and any i\, . . . , ik+m = 0,1, . . . ,d, 



(5.29) \\V n . . . V tk P t V lk+l . . . V lk+m fU < ct- d ^--^/ 2 \\f\\ 



■4,. 



Proof. We can apply [2H1 Corollary 2.17] to obtain that for each x £ M. , there 
exists a real-valued random variable ir x , depending on k and i\, . . . , ik+m, such 
that 

(5.30) V h ...V ik P t V ih+1 ...V ih+m f(x) =E[f(X?)ir x }. 
Furthermore, for each p £ [1, oo), there exists a constant C > with 

(5.31) sup E[\ir x \ p ] < c , t-deg(u,..,* fc+m )/2 > 

It follows that for p, q £ (1, oo) with 1 + i = 1, 
(5.32) 

\\v h . . . v lk p t v lk+l . . . v lk+m f\u, < sup ^r^i/P^)! • k x |] 

< \\fu sup ^xy^ixfyY'v ■n* x \ q ] 1/q 

x£R N 

< ct- deg(ii '-' ik + m)/2 \\f\y, 

where we apply Corollary 15.121 □ 



Corollary 5.14. Assume that the UFG and V0 conditions are satisfied. Then, for 
any mesh = t < ■ ■ ■ < t n = T and f £ C^(R N ), 

(5.33) \\PTf-Q tl -t ---Qt n -t n -JU 

< emu [(tn - ^-!) l/2 + e { \ T l :i )m/2 ) ■ 



CUBATURE METHODS FOR SPDES 



21 



Proof. We proceed as in the proof of [351 Proposition 3.6]. Assume first that / 6 
C|^ c (R Ar ); the general result then follows from a density argument. First, note that 

(5.34) 

\\PAtf ~ fh < SUp tKx^mfiXZt) - f(x)\] 

< sup ^(x)-^ sup \\Vf(sX* At + (1 - s)x)\\ ■ \\Xl t - x\\] 

x£f N s6[0,l] 

< ||V/|U sup ^(x)-^ sup i;(sXl t + (l-s)x) 2 ] 1 / 2 E[\\Xl t -x\\ 2 ] 1 / 2 . 

iei N sg[o,i] 

As ||sA^ t + (1 — s)x\\ < max(||A^ t ||, ||x||) for all s E [0, 1] and cosh is monotonic 
on [0, oo), we see that Corollary 15.121 yields 

(5-35) ||i'At/-/IU<C(A*) 1 /3||v/||^. 
By Theorem 15. 1 31 we obtain 

\\(PAt-QA t )p T -tfh< (a^)(™ +1 )/ 2 m v tl ...v lk p T ^ t f\u, 

(ii,...,i h )eAm 

(5.36) < C{At) (m+1)/2 (T - t)- m/2 \\Vf\\ 4 ,. 

Summing up in the usual manner, the claim follows. □ 

Corollary 5.15. Under the UFG and V0 assumptions, the cubature method con- 
verges of optimal order for f £ B^(M. N ) with V/ € B^(Wi N ) on graded meshes such 
as the ones chosen in [251 Example 3.7]. 

Proof. This follows directly from Corollarv l5 . 1 4l together with the density of Cg° (M. N ) 
in B^(R N ). D 

Appendix A. A counterexample 

Not every admissible weight function is D-admissible, as already the counter- 
example X = R, ip(x) :— 1 + x 2 + x~ 1 X(o,oo) w hh xa(x) :— 1 for x € A and for 
x (f: A the indicator of A from [T^l Remark 4.6] shows. However, such an assump- 
tion is necessary to be able to transfer differentiability properties to limits when 
using weighted supremum norms. 

Let us consider a concrete example. Choose the admissible weight function 
ip(x) := 1 + x 2 X(-oo.o) + x ~ 2 X(o,oo) on J = 1. Choose bounded, smooth functions 
f n , n e M, by 



(A.l) f n {x) 



0, x<0, 

1, x > n -1 , 



and monotone and smooth on (0, n _1 ) such that < f n (x) < 1 for all i£R. This 
can be done in such a way that \fn{x)\ < Cn on (0, n" 1 ) for some C > independent 
of n G N, for example by choosing fi as required and setting f n (x) '■— f\{nx). Then, 
for «eN and m > n, 

(A.2) ||/„ - / m ||. < 2 sup iPix)- 1 = 2(1 + n)- 2 and 

a;G(0,n- 1 ) 

(A.3) |/nki<2 sup i/>(x)- 1 Cn = 2C „ U 

from which 

(A.4) |/„-/mki<2C 



(l + n) 2 (1 + m) 2 
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It follows that (f n )neN is a Cauchy sequence in B{ (R). As evaluation functionals 
are continuous, we see that the only candidate for the limit is / = X(o,oo)- This 
function, however, is not in B^(R), and is not even continuous. 

This is not a contradiction to the characterization of (R) by continuity, as no 
set Kn < R} contains a neighborhood of x = 0. 
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